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Abstract 

From the open boundary t — J model, an impurity model is constructed in which 
magnetic impurities of arbitrary spins are coupled to the edges of the strongly 
correlated electron system. The boundary R matrices are given explicitly. The 
interaction parameters between magnetic impurities and electrons are related to 
the potentials of the impurities to preserve the integrability of the system. The 
Hamiltonian of the impurity model is diagonalized exactly. The integral equations of 
the ground state are derived and the ground state properties are discussed in details. 
We discuss also the string solutions of the Bethe ansatz equations, which describe the 
bound states of the charges and spins. By minimizing the thermodynamic potential 
we get the thermodynamic Bethe ansatz equations. The finite size correction of 
the free energy contributed by the magnetic impurities is obtained explicitly. The 
properties of the system at some special limits are discussed and the boundary 
bound states are obtained. 
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I. INTRODUCTION 



Strongly correlated systems have been of great interest over the past decade because of their 
importance as the most fundamental systems related to the theory of high- T c superconductors |l]J^] 
and these electron systems pose many theoretical challenges due to their essentially nonperturbative 
nature. The Hubbard chain and t — J model are frequently invoked as the basic models of the one 
dimensional strongly electron systems. When the repulsion energy U for two electrons located on 
the same atom is much larger than the bandwidth of the electrons, the Hubbard model reduces 
to the t — J model with an occupation of one electron per site. Then there are three states per 
site in the t — J model. Also the low-energy excitations of the three-band Hubbard model can 
be changed into an effective one-band t — J model || || . The t — J Hamiltonian is one of the 
most attractive model to explain high temperature superconductivity without any phonon pairing 
mechanism [^[0]. 

Recently great progress has been made for the one dimensional systems || , in particular for the 
quantum impurityproblems such as Kondo problem and tunneling in quantum wires. The theory of 
the Kondo effect || was developed in the early 1960s to explain the puzzle of the resistance of some 
metals, that is, the resistance start to increase as the metal is cooled below a certain temperature. 
This is due to that the local moments on the impurity atoms having an antiferromagnetic coupling 
to the spins of the conduction electrons and the coupling becomes stronger as the temperature falls. 
Now we know the impurities play an important role in the strongly correlated electron compounds. 
Even a small amount of defects may change the properties of the electron systems. 

Exacts results, even for a simplified model exploring just some features of the system, are always 
useful and provide a testing ground for approaches intended for the full problem of higher order 
of complexity. Therefore it is of interest to construct integrable systems of the strongly correlated 
electron systems including impurities. Indeed, there has been a long successful history to study the 
effects of the impurities in the many-body quantum systems within the framework of the integrable 
systems |o| Q, although the impurity usually destroys the integrability of the model when it is 
introduced to an exactly solved model. Andrei and Johannesson incorporated a magnetic impurity 
of the arbitrary spin into the isotropic spin-i Heisenberg chain with integrability preserved. The 
results were extended to the Babujian-Takhtajan spin chain in Refs. [^5|. Bedufig, E/3Zer and Frahm 
Jl6| ] solved the integrable model with the impurity coupled with periodic t — J chain Jl~7| — [3~0|] , in 
which the impurity is introduced through a local vertex. Schlottmann and Zvyagin introduced 
the impurity into supersymmetric t — J model [po[ via its scattering matrix with the itinerant 
electrons plp^ ] . The Hamiltonian of the system and other conserved currents can be constructed 
in principle by the transfer matrix. Zvyagin get that the low field impurity behavior for the periodic 
correlated electron chain coincides with that in the open chain up to mesoscopic corrections [2j| . 
They also discussed the magnetic impurities embedded in the Hubbard model |24|] and a finite 
concentration of magnetic impurities embedded in one-dimensional lattice via scattering matrices 
23 ]. The quantum impurity problem was also discussed by the use of the bosonization technique 

The study of completely integrable quantum spin open chains 1 28 - 32 1 is also an interesting 
subject. It deals with the system in a finite interval. The pioneering works of Cherednik and 
Sklyanin p^j28[ ] are the starting point for the one dimensional exactly solved models where there 
is a variant of the usual quantum inverse scattering method. The so-called reflection equations 
appeared as a new ingredient to the Yang-Baxter equations for describing the systems on a finite 
interval with independent boundary conditions on each end. Using this method, a lot of integrable 
models have arisen J3^-|3^] such as the Hofstadter problem and the reaction diffusion equations. 
The quantum-group-invariant transfer matrices and the Hamiltonians can be given by the certain 
limits of the open-boundary conditions p^- ^9| , ^ , ^l| ■ Some models have been solved with non- 
trivial generalizations of the methods developed by periodic and twisted boundary conditions. 

In 1993, Foerster and Karowski proposed the spl q (2, 1)— invariant t — J model with the 
quantum-group- invariant open boundary conditions and it was generalized to the SU q (n)- invariant 

chains by de Vega and Gonzalez-Ruiz in Ref. |4lJ . Gonzalez- Ruiz discussed also the open boundary 
supersymmetric t — J model commuting with the number operator n and S z [|42j . Recently, Zhou 
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and Batchelor [|43j studied the open boundary t — J model via an analytic treatment of the Bethe 
ansatz equations. Fan, Hou and Shi H] formulated the eigenvalues and eigenvectors of the t — J 
model with reflecting boundary conditions by the graded quantum inverse scattering method. The 
zero-temperature boundary effects in an open t — J chain with boundary fields were discussed by 
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Now we know that the quantum inverse scattering method provides us also a very useful tech- 
nique to construct the impurity model with the preserving of the integrability and this scheme 
has been adapted in Ref. ll^pl| , |2^ | . Very recently, the impurity model related to the t — J model 
was studied also in Ref. pqf But, when the impurities are embedded in the system with pe- 
riodic boundary conditions, some unphysical terms must be present in the Hamiltonian to keep 
the integrability though they may be irrelevant. Notice the fact that the impurities cut the one- 
dimensional system into 'pieces' when they are introduced and (then) the open boundary systems 

are formed with the impurities at the ends of the systems. In fact, Gonzalez- Ruiz |Q has pointed 
out that the boundary terms could take into account impurities or magnetic fields located at the 
boundaries in 1994. So, the integrable impurity models Jf7|-|52[ may be studied with the use of 
open boundary conditions. 

In this paper, we construct a Hamiltonian of the impurity model within the framework of the 
open boundary t — J chain. The magnetic impurities are coupled to the edges of the strongly 
correlated electron system and have arbitrary spins. The boundary R matrices are given explicitly 
and satisfy the reflection equation. They are compatible with the Yang-Baxter equation in the 
bulk. These guarantee the integrability of our model including the impurities with the arbitrary 
spins. The interaction parameters of the magnetic impurities with the electrons are parameterized 
by the constants which relate to the potentials of the impurities in order to contain the integrability 
of the system. The Hamiltonian with the impurities is diagonalized exactly by the use of the Bethe 
ansatz method and the Bethe ansatz equations are obtained explicitly via the quantum inverse 
scattering method. The integral equations of the ground state are derived and the properties of 
the ground state are discussed in detail. The finite size corrections of the ground state energy 
due to the impurities are obtained. The distribution functions of the rapidities are given. The 
string solutions of the Bethe ansatz equations are discussed and the coupled integral equations 
of the excited state are obtained. By minimizing the thermodynamical potential of the system 
with impurities, we get the thermodynamical Bethe ansatz equations. The finite size correction of 
the free energy is obtained for the two magnetic impurities. We study also the properties of the 
system in the high-temperature and the low temperature limits. The boundary bound states are 
constructed in our model. 

The paper is divided into eight sections. In section 2 we give the Hamiltonian of our impurity 
model where the generalized "permutation operator" is presented. In section 3 the integrability 
conditions are discussed and the boundary R matrices are obtained. Our Hamiltonian is diago- 
nalized exactly and the Bethe ansatz equations of the impurity model are obtained in section 4. 
In section 5 the coupled integral equations of the ground state are derived and the ground state 
properties of the system are discussed. In section 6, the string solutions of the Bethe ansatz equa- 
tions are discussed and the integral equations of the excited states are obtained. We get also the 
thermodynamic Bethe ansatz equations. The finite size correction of the free energy contributed 
by the two magnetic impurities is obtained. In section 7, the properties of the system at the 
high-temperature and the low temperature limits are discussed and the boundary bound states are 
constructed. Conclusions follow in section 8. In order to complete the four integrable cases of the 
t — J model, some boundary R matrices and the Bethe ansatz equations are given in Appendix A 
and B, respectively. 



II. THE HAMILTONIAN OF THE IMPURITY MODEL 

As is well known, the t—J model describes electrons with spin | on an one-dimensional lattice and 
permits the nearest- neighbor hopping (t) of the electrons. The large on-site Coulomb repulsion 
makes that the double-occupancy of every site impossible. There are two types of interactions 
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between the electrons on the nearest neighbor sites. One is the spin exchange interaction (J) and 
the other is the charge interaction (V ) independent of spin. Our starting point is the following 
Hamiltonian 

G-1 G-1 G-1 

ff = -«EE (c£ c j+la + c+ +lCT c oa ) + j E s 3 • + v E 

3=1 £T=TJ. 3=1 3 = 1 

+J a Si • d a + V a ni + J b S G ■ d b + V b n G , (1) 

where C^.(Cj CT ) is the creation (annihilation) operator of a conduction electron with spin a on 
the site j; Sj = | ct , C^ a aM iCj a i is the spin operator of the conduction electron; d a ^ are the 
spin operators of the impurities; J a . b , V ab are the interaction constants of the impurities with the 
electrons and the scattering potentials of the magnetic impurities , respectively; rij — C^Cj-\ + 
CjiCjiis the number operator of the conduction electrons; G is the length ( or site number ) of 
the system. For the bulk of the system, the scattering matrices S of the electrons satisfy the 
Yang-Baxter equation: 

Si 2 (ki, fa)Sx 3 (ki, fa)S 23 (fa, fa) = S 23 (fa, fa)Si 3 (k 1 , fa)Si 2 (fa, fa) (2) 

which gives the factorization condition for the integrability of the free boundar y a nd periodic 
boundary t — J models. For the open boundary system, the reflection equation |3^ , |28[| should also 
be satisfied, i.e. 

Si 2 (fa, fa) R (fa)S 12 (fa, ~fa) R {fa) =R (fa)S 12 (fa,-fa) R (fa)S 12 (fa, fa), (3) 
where the boundary R matrices are written down as 

1 2 

R(fa) = R(fa) (g) id V2 , R (fa) = idvi <8> R(fa) 

for matrix R 6 End(V). 

To consider magnetic impurities of arbitrary spins, we define an operator as 

Pd = = = = (1 + 45 • d) (4) 

y/l + 41(1 + 1) 

where S and d are the spin operators of the electron and impurity, respectively. It has the property 
P% — 1 and degenerates to the ordinary permutation operator when the impurity spin d takes its 
value I = 1/2. We call Pd the generalized "permutation operator". The Hamiltonian may be 
written in the form 

N N 

# = - E ( T , + + T r) + E fa + + K bj s XjtG ) (5) 



V 3 3 / 1 Z-^i 

3=1 3=1 



where we have set t = 1. The translation operators Tf are defined as fc3j 
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T j • • • ,xj, ■ ■ ■ ,x N ) = $>(xi, ■■■,Xj±l,'- ■ ,x N ) 

where ^(x±, ■ • • ,Xj, ■ ■ ■ , £/v) is the wave function of N conduction electrons. The operator Kj acts 
on the wave function ^(xi, ■ • • , Xj, • • • , xn) in the form 

N 

Kjfy(xi,x 2 , ■ ■ ■ ,x N ) = >^ 6 Xj>Xt+ tKij9(xi, x 2 , ■ ■ ■ ,x N ) 

i=l 

where Kij — V — j + j[Pij denotes the interactions between the conduction electrons. The 
permutation operator Py permutes the spins of the i -th and j-th electrons. The interactions 
between the impurities with arbitrary spins and the electrons are given by the relations 
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K aj = V a f + 



J a J aV /l + 4s a (s a + l) 



Abj = Vfe - — + Aj 



with 



v /l + 4s a (s a + l) 



(1 + iSj ■ da) 



(6) 



Pbj = 



^1 + 4^6(36 + 1) 



(1 + 45, • 4) 



(7) 



where Sj are the spin operators of the conduction electrons; d a and db are the spin operators of 
the magnetic impurities with the arbitrary spin values s a and s&, respectively. The boundary R 
matrices of our impurity model will be given explicitly in the following section. 



III. THE BOUNDARY R MATRIX AND INTEGRABILITY CONDITIONS 

The boundary R matrices satisfy the reflecting Yang-Baxter equation for the open boundary 
system and the reflection equation is compatible with the factorization condition. These give the 
integrability of the one- dimensional many body system. In our model, the boundary R matrix 
should take the form 



R = cxp(if) 



q-ic-i(2l + l) P d /2 
q + lC + l(2l + l)P d /2 , 



(8) 



where Pd is the generalized "permutation operator" , q = | tan | and c is an arbitrary constant. 
Here we have set J = 2 and V = 3/2. The other integrable cases will be given in the Appendix. 
Put K a (tyj = m + nPd ■ From the boundary Yang-Baxter relation we have that 



2n 



21 + 1 



21 + 1 



2l + l\ k 

— ^ — tan x 
2) 2 



Q 



(to — l) 2 — n 2 



tan 



(to + 1) 2 -n 2 



= 0. 



(9) 



This relation gives the restricting condition imposed on the coupling constants of our impurity 
model. It should be satisfied to maintain the integrability of the system when the magnetic 
impurities with the arbitrary spins are introduced to the t — J model by coupling the impurities 
to the edges of the system. From the above relation, we have that 



rn = 



(2l + l) z 



4(/±c)(/ + c+l)' 



2Z + 1 



2(Z±c)(Z + c+l)' 



The phase factor in the boundary R matrix related the magnetic impurity at the left edge of the 
system is given by 
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r . (]A] _ 2 J a (q + iCq) expjik) + i[A + (4V a - J a ) exp(ik)} 
exp L t VoW J - 2Jq((? _ iCa) eM _ tk) _ . [4 + (4K _ Jq) exp( _ ifc)] • UUj 

The phase factor exp [iip (k)\, corresponding to the scattering of the electron with the impurity at 
the right edge of the system, is given by the substituting b for a. Without loss of generality, the 
interaction parameters J a ,b and V a ,b take the following forms: 

2 

(s a + c ) (s Q - c a + 1)' 11 

y - ^-(2 5a + if -3 

4( Sa + c a )( So - Co + l) 1 j 

The parameters J b and are given similarly by the substitution of index a by index b in the 
above relations. By the use of these parameterizations of the interaction constants J a ,b and V Qjb , 
the phase factors take the form: exp [iip a (k)] — exp [iipj, (k)] = — exp (ik). Then, the boundary R 
matrices for the impurity model with the arbitrary spins can be expressed as 

, s 2g 7 - - 2ic a - i (2s a + 1) P ai , s 

«.(*,,,) (13) 

of i \ r •; /o^ , im 2( ?j ^ 2ic b - i (2s b + l)P bj 

Rb(-kj,(Tj) = - exp -ifc ■ 2G + 1) — ,n P 14 

2gj + 2iCb + i (2sb + 1) Pbj 

with <7j = i tan , respectively. Here the operators P a j and Pbj are the generalized "permutation 
operators" and have the expressions (g) and (|7]). 

IV. BETHE ANSATZ EQUATIONS OF THE IMPURITY MODEL 

Firstly, we can write down the wave function 9 ait<Tai ... „ K (xi, X2, • ■ ■ , Xjy) in the region < Xq\ < 
XQ2 < ■ ■ < xqn < G — 1 as 

^a 1 ,a 2 ,---,<J N {xi,X2, • ■ ■ , Ejv) 

= X] X] ep£rA CT<31 , CTQ3 ,..., CTQN (rpQifepQi,rpQ 2 fepQ2,---,'"FQivfepQiv) 

JV 

•exp^^rpjfcpj-ar^], (15) 
3=1 

where the coefficients A tTQ1 , aQ2 ^..., aQN (rpQikpQi,rpQ 2 kpQ2, - ■ ■ ,r P Q N k P Q N ) are also dependent 
on the spins of magnetic impurities which are suppressed for brevityand sp — 1(— l),when the 
parity of P is even(odd), e r — Yij=i r m which r takes the value+lor — 1. The coefficients A are 
related by the scattering matrices. The boundary R matrix satisfies the reflecting Yang-Baxter 
equation (^ ). The scattering matrices in the bulk satisfy the Yang-Baxter equation (|| ). By 
the use of the standard Bethe ansatz method, we can diagonalize the Hamiltonian with the two 
magnetic impurities |ll|Q. The boundary R matrices can be cast into the forms: 



| gj - icq + i jsg + 1/2) S jo {kj,ko) 
qj + ic a - i (s a + 1/2) Sj (-kj,k ) ' 



R a ( kj ^) = - exp(^) ? -J^Ll ■ ; tL^ p^222 ^ , (16) 



or u ^ r -u ion i 1 M * _ ic& + * ( Sfe + 1 / 2 ) S 3N+i( k 3> k N+i) n „, 

Rbi-k^CJj) = - exp -zfc ? - (2G + 1) — : — 77777 t — 7 ; r, (17) 

3 J 3 q 3 +ic b - i(s b + 1/2) S jN+1 (-k j ,k N+i y 



G 



where 

sAkM = ^ - <io -Hs a + 1/2) p 

3 V J 1 qj -q + i(s a + 1/2) V ' 

c tu j. \ q 3 -qN+i-i(s b + 1/2) P b3 

b jN +i{Kj,k N+ i) —— , , (19) 

?j - qN+i + i {sb + 1/2) 

with go = «c a and qN+i = icb- Notice that the scattering matrix in the bulk has the form 

S jl (k j ,k l )=- qj ~ qi ~f j \ = 1,2,- --,N) (20) 
qj-qi+i 

where qj — | tan%. The present case ( J = 2 and = 3/2 ) corresponds to triplet scattering and 
not to the ordinary t — J model. The symmetry of this model is 51/(3) , rather than a graded FFB 
super-algebra which corresponds to the traditional i—J model. In this case |l7|], x* = V/2+J/8 = 1 
and x s = V/2 — 3 J/8 = 0. These quantities are interchanged in Ref. ]l7j ], but this was corrected 
in Ref. g. If we set 

5 Ti (A) = - A ~ g '~5 T ' : (i = 1,2,.--,JN0 
A - qi + i 

a m 2X-2ic a -i(2s a + l)P ar 

Or0(AJ = — — — — jz — — t , (21) 

2A - 2ic a + i (2s a + 1) 



2A-2^c fc -i(2s b + l)P fcr 
2A - 2ic b + i (2s b + 1) ' 



S t n+i (A) 
then we get the equation 

Tr \T (A) T- 1 (-A)l L 4> = exp {2ik,G) \ - gj gj + ' Ca ~ \ + ^ 
L 1 ^ V yJlA= ^ V J ' i/2 + q j q j -ic a +i(s a + 1/2) 

Qj +ic b -i{s b + 1/2) ^ . 

qj - ic b + i (s b + 1/2) 

where $ is the eigenstate of the system and T(A) is defined as 

T (A) = S TJ (A) S T0 (A) 5 Tl (A) • ■ ■ 5 TJ _! (A) 5 TJ+1 (A) • ■ ■ 5 tA t +1 (A) . (23) 

Then we have the following Bethe ansatz equations: 

M 

e (2^) 2G Y[ e (2qj - 2X?) e {2q 3 + 2A„) 

2q 3 \ ( 2q 3 
2c a + 2s a + lJ \2c b + 2s b + l 

N 

■ II e fe-?0efe+9i)> (24) 

Ar» \ / X n \ ( X 



s Sa + Co / V s fc + °b J \S a — C a J \S b ~ Cb , 

N 

■\[e(2X a -2qi)e{2X a + 2qi) 

M 

= e (Aa-A^)e(A a + A /3 ), (25) 
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where we have used the notation e(x) = (x + i) / (x — i) J55|. Notice that e(±oo) = 1. For the case 
of J = —2 and V — —3/2, only the energy changes sign. The Bethe ansatz equations for this case 
are the same as for antiferromagnetic coupling. 



V. GROUND STATE PROPERTIES 



The eigenvalue of the Hamiltonian can be written down as 



E = 2N - 



N 

£ 

3 



q 2 

1 H 3 



(26) 



1 k ■ 

with qj — g tan . In order to obtain the ground state properties for the case of the triplet 
state scattering we should find out the solutions of the above Bethe ansatz equations. The sets 
of the rapidities {qj} and {Aq,} have real and the complex solutions. Complex solutions of the A a 
correspond to the excited states and will be discussed in section 6. 

For J — 2 and V = §, the model is the SU(3) invariant t — J model where electrons in triplet 
states are scattered but not those in singlet states. For the bulk this model is isomorphic to the 
spin 1 Heisenberg chain with SU(3) invariance. It has no graded super-algebra and corresponds 
to three bosonic degrees of freedom. All rapidities are real for the ground state. Now we take the 
thermodynamic limit and introduce the distribution functions p (q) and a (A) corresponding to the 
parameters q and A. The corresponding ones of the holes are denoted by p h (q) and a h (A). So the 
integral equations describing the ground state of the impurity model can be written as 



a { q, l) + \ J dX<7 ^ 



" ' q ~ A ' 2 7 +a \ q + A ' \ 



p(q) + p h (q) + ^ 



a ( q, c a + s a + i j + a (q, c b + s b + ~ 



a\q,~ 



+ 2 / P W) d l' i a (l - l\ l)+a(q + q', 1) 



(27) 



1 

2G 



a (A, s a + c a ) + a (A, Sb + c b ) + a (A, s a - c a ) + a (A, s& - c b ) + a ( A, - 



-\ J dqp [q) 



'' ' ' X ~ q ' 5) +a \ X + q, \ 



1 



a (A) + a h (A) + - / dX'a (A') [a (A - A', 1) + a (A + A', 1)] 



(28) 



where a (A, rf) = ^1 2 - The ground state properties can be obtained from the above coupled 
integral equations. The details are as follows. 

By taking into account of the above distribution functions, the energy per site of the system can 
be expressed as 



G=^- 2n J dqp{q)a { q >2) 
for the case of J = 2 and V = 3/2. The number of the particles per site is 

N f 

■q = dqp (q) . 
The magnetization of the impurity model can be written as 



(29) 



(30) 
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^ = \Jd qP {q)- jdXa(X) + 1 -^ (31) 

where l a = 1— s a , 2— s a , ■ • • , s a and 4 = 1 — s&, 2— ■ • • , s&. s a and Sb are the spins of the impurities 
located at the two ends of the one-dimensional lattice system. The Fourier transformations of the 
Bethe ansatz equations (^7|) and (|2^ ) give that 

e - M/2 + e -M/ 2 ~ (w) = ^ + e - H ) + ^ (w) + _L ?G (w) f (32 ) 

P M e-l^l/ 2 + ^ G (w) - 5 (w) (l + e-l-l) + a" M , (33) 

where the tilde denotes the Fourier transform and the functions pa (w) and ctg (w) are the Fourier 
transformations of p G (g) and oq (A) given by 

PG («?) = a [ ?,c a +s a + - ) + a ( g,c b + s b + - ) - a (g, - ] , (34) 



aa (A) = a (A, s a + c a ) + a (A, s b + c b ) + a (A, s a - c a ) 

+a(A, s& - Cb) + a ( A, - ) , (35) 



which are the finite size corrections due to the magnetic impurities and the open boundary con- 
ditions. Now we first consider the half-filled band of the ferromagnetic case of a (A) = 0. The 
nonzero distribution functions are 

pH= 2c^hf -2GTT^M' (36) 

and the energy is 

E/G = l-2ln2. (38) 
The finite size correction to the energy in this case is 

contributed by the impurities, and the open boundary term gives the correction — In 2. For the 
nonmagnetic case, when the number per site is 2/3, the distribution functions take the form: 

2 cosh- 1 e M/ 2 r 
= TT2^L + 2GU2^ \° G < W > - 2 ~ PG ^ C ° Sh 2 I ' (40) 



_ li e l"l /2 r ~ w i 

cr (w) = ! < 2er G (uj) cosh p G (uj) > (41) 

V 7 l + 2coshw 2Gl + 2coshw I GV 7 2 PG V '! K ' 



1 1 e N/ 2 

l + 2coshw 2Gl + 2coshw 
with the energy being 

E 4 tt\/3 



In 3, (42) 



G 3 9 
and the finite size correction of the energy is 

2 7 rv / 3 1 f°° a G {uj) - 2p G (uj) cosh f 
~9 2 i_ 00 1 + 2 cosh w 

This is the contribution of the magnetic impurities and the open boundary terms give the correction 
-2ttV3/9. 



, Z7TV o 1 / o G (UJ - Lp G w tusu 2", /.ON 
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VI. EXCITED STATE 



The eigenstates of the model are specified by two sets of rapidities, {qj} and {Xj}, for the charges 
and the spins, respectively. In this section, we will also restrict the discussions on the case of J = 2 
and V — 3/2. In the thermodynamic limit these rapidities are classified according to the string 
hypothesis: |22] , |56|] 

Q2' j =02 + £(n + l-2j), 3 = l,V",n 



A^' = A^ + -(n+l-2j), j = l,2,.-.,n (44) 

where n — 1,2, ■••,+oo. These strings describe charge and spin boundstates, respectively. The 
real parameters g™ and Xp are related to the momentum of the center of mass of the boundstate. 
Their distribution functions are p n (q), a n (A) and the ones corresponding to the holes are p 1 ^ (q) 
and cr^ (A) . By using the notations 

[n]f(k)= a(k-k',^jf(k')dk' forn^O, (45) 
[0}f(k) = f(k), 

the integral equations can be cast into the following compact forms: 

co m-1 _^ 
fl H) + EE^™ +1 + 2 3\ °™ («) + («) 

m=l 3=0 



= Pn (?) + 51 4»mA" («) ' ( 46 ) 
m— 1 

- oc m — 1 OO 

— o£ (A) + £ ^ [n - m + 1 + 2j] Pm (A) = a h n (A) + £ A„ m a m (A) (47) 

m—l j — m— 1 

where the terms with the factor 1 / (2G) are the contributions provided by the magnetic impurities 
and the boundary terms. The operator A nm is defined as 

A nm = [\n - m\] + 2 [\n - m\ + 2] + 2 [\n - m\ + 4] H h 2 [n + m - 2] + [n + to] (48) 

and 

n-l 

P« (?) = a (q, |J - 5] {« (?) «a + c a - j + |) + a (g, s 6 + c b - 3 + ^) j , (49) 
3=0 

n\ v—v /, n — 1 \ /, n — 1 



CT n ( A ) = a ( A ' 2 ) + 5Z 1 a ( A ' S °-+ C a- j ^ — ) + a ( A, S b + C b - j + 



+a ( A, s a - c a - j + n ] + a ( A, s h - c b - j + - 1 I ^ (50) 
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Now the number of the electrons per site is expressed by 

q = E m I p™ dq - ( 51 ) 

m— 1 

By minimizing the thermodynamic potential we get the thermodynamical Bethe ansatz equations: 

n (2 - A - H) - 2na (k, § ) , , ^ 

1 - In (1 + + E nm ln C 1 + C, 1 ) 

m— 1 

00 n — 1 

- E E t m - n + 1 + ^ ln C 1 + r >™) = °- ( 52 ) 

m— 1 j— 



2i?n 



ln(l + 77„)+ E A 

It m 

In (l + V7 m x ) 



rn— 1 



where we have set that 



OO 71 — 1 

- E E[ m - n+1 + 2 ^ ln ( 1 +^ 1 ) = ' ( 53 ) 

m— 1 J=Q 

Of course, the above expressions can be changed into the several other forms as the ordinary system 
in one dimensional lattice system. The impurities located at the ends of the system do not change 
the free energy for the bulk |22|]5rf l . The finite size correction of the free energy contributed by the 
boundary term is 

F bou = ~ £ W In {1 + C 1 (0)} - § E N ln i 1 + Vn 1 (0)} ■ (54) 

n— 1 n— 1 

The free energy due to the two magnetic impurities with the spins s a ,b is 

t 00 r 

F imp = ~ 2 E J d ^Pn,im P ^ (l + C') - 2# (l a + h) 



71=1 
OO 



- E J d\a° imp ln (1 + 77,; 1 ) , (55) 



n=l 

where l a = 1 - s a , 2 — s a , ■ ■ ■ , s a ; k = 1 - s&, 2 - s fc , • • • , s fc ; and 

n-l 



1)1.1 



M = - E { a Sa + Ca ~ J + f ) + a ( 9 ' Sb + Cb ~ i : + \ ) } ' ( 56 ) 



3=0 



71—1/ -, 

71—1 



<im P (A)= E E a ( A ' Si+rC '--?' + ^-) • ( 5? ) 
r=±l,J=a,6 j=0 ^ ' 

The interesting observation that the finite size correction due to the magnetic impurities is finite 
even if the temperature T — ► 0. The properties of the system at the high-temperature limit and 
the low temperature limit can be discussed by the use of the thermal Bethe ansatz equations. And 
the boundary bound states exist in our model. These will be studied in the following section. 
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VII. SPECIAL LIMITS AND THE BOUNDARY BOUND STATES 



In order to discuss the low temperature properties of the model and the impurity free energy at 
high-temperature, we may change the thermal Bethe ansatz equations into the forms: 

2tt ([01 + \2Vj~ 1 a(k i) ~ 

U J [ 1! ^^+Gln-±-f-+ln6=0, (58) 

1 1 + ?2 



- U „-l l+ „ + l 

ln£ n = Gln— — — r — — — r -, 77 >1, 

(! + ? 7n-l)( 1 + ?7n+l) 



(59) 



2^([0] + [2])- 1 a(fc,l) ~ 

G In (^2^2) + In (6m) = 0, 



(1 + C-1) (1 + C+1) 



(60) 
(61) 



lim [n + 1] In 



1 + 6 



n] In 



1 + 6+i \ A-2 + H/2 



(62) 



lim I [n + 11 In 1 + \ - [77] In 1 + ^"t 1 ) = -— (63) 




with the use of the formulae A iOT - GA 2m = 5 lm , A nm - G (j4„_i m 4- A n+1 m ) = 5 nm [n > 1), and 
[n + 1] A nm - [n] A n+X m = if m < n or - [m + 1] - [771 - 1] if m > n where G = [l]/([0] + [2]). 
Furthermore, we have that linin^oo (In 6) / 71 = (2 — A — H/2)/T and lim n ^oo (In77„) /n = if/T 
where H is the magnetic field and A is the chemical potential of the system. Define the symbol 
function sign(a) as 

sign(a) = 

Then, in the high-temperature limit, the free energy of the impurities is 
F imp = — In 2 Lign(s l + c l + ~+r 1 ) 

- ^ sign(si +r 2 ci + n) > . (64) 

r 2 =±l J 

The boundary term give the contribution f& ou = — 3Tln2. The above result means that the entropy 
of the model related to the impurities is dependent on the interacting parameters of the impurities 
with the electrons and the spins of the impurities. At the low temperature limit, from the thermal 



Bethe ansatz equations, we get thatJ57|]58| the functions 6 have the different asymptotic properties 
for 77 < 2 and 77 > 2. It means that |57|] / = 2 and this gives that r = 4/ (/ + 2) = 1. By considering 
the free energy expression in the bulk, we find that at low temperature, the specific heat behaves 
as Cbui ~ T. Since the asymptotic property of the term ln(l + 77" 1 ) are different from the one of 
the term ln(l + 67 ) & t the low temperature limit, the impurities and the open boundary term 
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contributions to the specific heat at low temperature are nonlinear at the temperature. Although 
it is difficult to get the exact expression of the specific heat, we know that this nonlinear relation 
with the temperature is related to the interacting parameters and the spins of the impurities. 

The boundary bound states can be formed in the present model in the charge sectors. When 
\s a + c a \ < i, the boundary bound state can be denoted by 

q 3 ' = ±i (s a + c a + ~ 

q^=^Us a + c a -\). ( (ir, i 

3.2 if. 3 



q ~F2\~ u, ~ Lt '2 
in the charge sectors. It carries the energy: 

(3) = 3(2s a + 2c a + l) 2 {2s 2 a + 2c\ + As a c a + 2s a + 2c a - 5) 

ed9 (s + c a )(s a + c a + l)(2s + 2c a -3)(2 Sa + 2c + 5) ' 1 ' 

Obviously, when we substitute the subindex a by b, we also get the bound state with three imaginary 
modes localized at the end of the system. The total momentum of the bound states are zero. 
Another kind of the boundary bound state with five imaginary modes can be expressed as 

q 5 ' = ±* (s a + c a + \ 

9 5 ' 1 = (3s a + 3c a + 1 

q 5a = Tt; (3s a + 3c Q + | ) , (67) 



q 5 ' 3 = ±i \ .s 
q 5 ' 4 = ±i 



2 V 2 
3 
2 
1 
2 



when — | < s a + c a < — i . It has the energy 



n(5) = 5 (72{s a + c a ) 4 + 144(s a + c a f - 58(.s a + c a ) 2 - 130(g a + c a ) + ll) 
' ed9_ (s a + c a -l)(s Q + c Q + 2)(6s a + 6c a -l)(6sa + 6c a + 7) 



The imaginary modes with the subindex b also form the boundary bound state. The momentum 
of the state, which is the sum of the ones of the imaginary modes, is also zero. These bound states 
are at the ends of the system. 



VIII. CONCLUDING REMARKS 

Recently, in an interesting work J2j| , Zvyagin get that the low-energy magnetic behaviors of an 
impurity in a chain with periodic boundary conditions and with open boundary conditions coincide 
up to mesoscopic corrections of order of G _1 and it was shown that this property is independent of 
the impurity position in the open chain. From the above discussions, we get the same conclusion. 
The impurity part of the Hamiltonian for the open boundary condition in Ref. |p3| is 

H lmp = [K aP /(l + ul)](J% jf + H.c.) (69) 
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with the "weak link" (1 + Uq) 1 < 1 where K a p = StrJ a J@ and Jq± '* "' are the generators of 
the supersymmetric algebra sZ(l|2). In our model, for the open boundary conditions, the impurity 
part of the Hamiltonian is 

H imp = JaSi ■ d a + V a ni + J b S G • d fc + V b n G . (70) 

J a .b and V a b depend on the interacting parameters c a .t and the spins s a .b of the impurities (See 
relations jl^ - 12 )). The finite size corrections of the ground state due to the impurities are 



obtained and the properties of the system at the high-temperature limit and the low temperature 
limit are studied in this paper. We also discussed the boundary bound states which are localized 
at the edges of the system. 

To conclude, we have constructed from the open boundary t — J model a Hamiltonian with the 
two magnetic impurities which have arbitrary spins . The impurities are coupled to the ends of the 
system. The interaction parameters of the impurities with the electrons are parameterized by the 
constants which related on the potentials of the impurities in order to preserve the integrability of 
the system. The boundary R matrices of the impurity model are given explicitly and satisfy the 
reflection equation. In fact, our impurity model has also the following R matrices, 



R a (k,,a,) = - exp [iki T i0 a (tt - kj)] ^ ^ - f Sa + \ ) Paj , (71) 

1 ° V JH 2 qj +2i Ca + l(2Sa + l)Pa 3 



R b (-k 3 ,a 3 ) = - exp [-i kj (2G + 1) T Uh, (* - kj)] ^' f\ j f' b ± ]\ g 3 (72) 

Zqj + iicb + i \2,Sb + 1 J -t%j 



with qj = | tan -j- . Here the interaction parameters are denoted by 

j a = ? , 

(Sa ~ C a ) (s a + C a + 1) ' 



4 (s a - C a ) (s a + C a + 1) 

The function 9 is defined as 

( _ 1 (2c^ - 2sj - 2s a ) cos k + 2sj - 2c\ + 2s a + 1 + 2ic a sin k 
a ^ ' ~~ i 11 (2 C 2 - 2s 2 a - 2s a ) cos k + 2s\ - 2c\ + 2s a + 1 - 2ic a sin k ' 

The parameters Jb, Vb and the function 9b (k) can be obtained from the above relations with the 
substitutions of the index b for a . By comparing the expressions of J a ,b and V a ,b with the ones in 
section 3, we know that the different expressions of the boundary R matrices are because of the 
transformations c a — > — c a and — > — c^. This reminds us the R matrices of the impurity model 
can be written down as the other forms with similar substitutions. 

The Hamiltonian including the arbitrary spin impurities is diagonalized exactly by the Bethe 
ansatz method. The coupled integral equations of the ground state are derived and the ground 
state properties are discussed in details. The distribution functions of the rapidities are given 
also. Furthermore, we study the string solution of the Bethe ansatz equations. The integral 
equations of the excited state are obtained. By minimizing the thermodynamic potential, we 
get the thermodynamic Bethe ansatz equations. The finite size correction due to the magnetic 
impurities is obtained for the free energy. It is interesting that this correction is finite when the 
temperature T — > 0, but have the value related to the magnetization of the impurities and the 
external magnetic field. We discussed also the properties of the system at the high-temperature 
limit and the low temperature limit. Two kinds of the boundary bound states are constructed 
also for the case of J = 2 and V — 3/2. Finally, we point out the impurity contributions are 
derived in the above discussion under the thermodynamic limits since the Bethe ansatz equations 
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can be changed into the coupled integral equations and there solved by a Fourier transform. So, 
it is an interesting problem to discuss the impurity effects for the finite lattice case, de Vega and 
Woynarovich fl59|] have given a method to calculate the leading-order finite size corrections to the 



ground state energy. By using the finite size scaling technique 60 - 62 1 , the case of some other 
integrable models such as Hubbard model and Heisenberg spin chain are also studied |]63|-|65|| . It is 
worth studying the impurity effects in different sectors of the present model further and its critical 
properties are also interesting topics for further discussions. 



APPENDIX A. THE BOUNDARY R MATRICES 

When the interacting parameters J a and V a j, take the forms 

2 

Ja = (s a - Ca )( Sa + c a + iy {AA) 

4c^-(2 So + l) 2 + l 
Va A{s a ~c a ){ Sa + c a + l) 1 ' ! 

for J = 2, V — —1/2, the boundary R matrices for the impurity model with the arbitrary spins 
can be written down as 

(A - 3) 



with q = | cot | . When J = — 2 and V = 1/2, J a and V a are 

T = - 

[s a ~ Co) (S a + C a + 1) 



Ja = —< TT—> ~TT' ( A - 5 ) 



V a = - A f a ~ ^ 2S ° + / + 1 (A.6) 



4(s a - c a ) (s a + c Q + 1) 
When J = — 2 and = —3/2, we have that 

^° = — 7 I w TTT' (A-7) 

= 4cg-(2 S q + l) 2 -3 
a 4( So + c„)( S o-c + l)- 1 j 

The boundary R matrices are same as the relations (|l^ - |l4|) . 

APPENDIX B. THE BETHE ANSATZ EQUATIONS 

From the standard Bethe ansatz procedure, we get the following Bethe ansatz equations: 
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/ot, \ COt if + 2iC " + * ( 2s « + !) C0t ^ + 2iC & + * ( 2g b + 1) 

exp \2Gikj) # r: 

cot -j- — 2ic a — i (2s a + 1) cot -4- — 2zc& — i (2s& + 1) 



-| — r COt 

n 



2A, 



i cot 



2A, 



=i cot ^ 



2A, 



% cot -j- + 2A/3 



(B.l) 



A a + zc a + is a X a + icb + ist A 



icb + isb 



X a + ic a — is a X a + icb — isb X a — ic a — is a X a — icb — isb 



N 

n 

i=i 



2X a -cotf+i 2X a + cot f + i 



2X C 



cot • 



i 2X r 



cot • 



M . . . . . 

= II A a "x g+1 A a t^ + 1 («=1.2,-,W) (B.2) 
A Q — A/3 — z A Q + Xg — l 

/3=l(/3#a) 

for the case of J = 2 and V = —1/2. When J = — 2 and V = 1/2, the Bethe ansatz equations are 
the same as the above ones, but the energy ch ang es sign. The Bethe ansatz equations for the case 
of J = — 2 and V = —3/2 are the same as (|2^-p5|) and the energy changes sign also. 
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